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Abstract. We discuss the first astronomical detection of the CF^ (fluoromethylidynium) ion, obtained by observations of the 
J = l-O (102.6 GHz), y = 2 - 1 (205.2 GHz) and 7 = 3 - 2 (307.7 GHz) rotational transitions toward the Orion Bar region. 
Our search for CF+ - carried out using the IRAM 30m and APEX 12m telescopes ~ was motivated by recent theoretical models 
that predict CF^ abundances of few x 10"'" in UV-irradiated molecular regions where C^ is present. The CF^ ion is produced 
by exothermic reactions of C^ with HF. Because fluorine atoms can react exothermically with H2, HF is predicted to be the 
dominant reservoir of fluorine, not only in well-shielded regions but also in the surface layers of molecular clouds where the 
C^ abundance is large. The observed CF^ line intensities imply the presence of CF^ column densities > 10'^ cm"^ over a 
region of size i 1', in good agreement with theoretical predictions. They provide support for our current theories of interstellar 
fluorine chemistry, which suggest that hydrogen fluoride should be ubiquitous in interstellar gas clouds and widely detectable 
in absorption by future satellite and airborne observatories. 
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1. Introduction 

Of the ~ 120 interstellar molecules that have been detected 
to date' (e.g. Miiller et al. 2005), more than 95% involve just 
six elements - hydrogen, carbon, nitrogen, oxygen, silicon, and 
sulfur - all of which have moderately high abundances (> 10"^ 
relative to hydrogen in the Sun). Of the molecules that are 
formed in the interstellar medium, only four species contain- 
ing other elements have been reported previously: HF (Neufeld 
et al. 1997), FN (Turner & Bally 1987), HCl (Blake, Keene & 
Phillips 1987) and FeO (Walmsley et al. 2002). In the case of 
those elements of relatively low abundance (F, CI, and P with 
solar abundances of 3 x 10"^, 3 x 10"^ and 3 x 10"^ respec- 
tively with respect to H), the detected molecules are all stable 
species with large dissociation energies (D^ = 5.87, 4.43 and 
7.05 eV respectively for HF, HCl and PN^). Not surprisingly, 
our understanding of the chemistries of interstellar fluorine-, 



' In addition, roughly two dozen molecules have been detected 
exclusively in circumstellar gas - specifically around the source 
IRC -I- 10216. These molecules include species containing the ele- 
ments Na, Mg, K, and Al, in addition to the elements represented by 
molecules formed in the interstellar medium. 

^ Derived from thermochemical data recommended in NIST 
Standard Reference Database Number 69 (a.k.a. the NIST Chemistry 
Web book), Eds. P J. Linstrom and W. G. Mallard, March 2003, 
National Institute of Standards and Technology, Gaithersburg MD, 
20899; available on-line at http://webbook.nist.gov_ 



chlorine-, phosphorus- and iron-bearing molecules is limited 
by the paucity of detected species that might help constrain 
theoretical models. 



Notwithstanding this lack of observational constraints, the 
chemistry of fluorine-bearing molecules in the interstellar 
medium has been the subject of a recent theoretical study 
(Neufeld, Wolfire & Schilke 2005a; hereafter NWS). The prin- 
cipal conclusions of that investigation were: 1) that hydro- 
gen fluoride is formed rapidly by the exothermic reaction of 
F atoms with H2 and becomes the dominant reservoir of gas- 
phase F nuclei over a wide range of conditions; (2) that HF is 
abundant even close to UV-uradiated cloud surfaces where C^ 
is the dominant reservoir of gas-phase carbon; and (3) that reac- 
tion of HF with C^ can lead to potentially-measurable column 
densities of CF^, a very stable molecule (dissociation energy, 
Dq = 7.71 eV) that is isoelectronic with CO. 

NWS showed that the chemistry of interstellar fluorine is 
qualitatively different from that of other elements because flu- 
orine atoms can react exothermically with H2, the dissociation 
energy of hydrogen fluoride being the largest of any neutral 
diatomic hydride and HF being the only such molecule more 
strongly bound than H2. Fluorine is therefore unique among 
the elements in having a neutral atom that reacts exothermi- 
cally with H2. 
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Fig. 1. CF"^ 7=1-0 spectra observed at each of the positions in 
the Orion Bar, with the observed position for each observation 
indicated on a '^CO 7 = 3-2 map of the source (grayscale map, 
from Lis & Schilke 2003). For the strip centered on (a, 6) - 
(05h 35m 22.8s, -05" 25' 01") (J2000) (the "Orion Bar (CO)" 
position), the integrated CF+ J = 1-0 antenna temperatures are 
33±6, 45+8, 86±10, 125+13 and 57±8 mKkm s'^ from NW to 
SE (quoted uncertainties are 1 tr statistical errors); for the strip 
centered on (05h 35m 25.3s, -05" 24' 34") (J2000) (the "Orion 
Bar (HCN)" position), the corresponding values are < 19.5, 
< 15.8, 72 ±8, 119± 15,andl40±30mKkms-i. 

Motivated by the predictions of NWS, we used the 30m 
TRAM telescope and the 12m APEX^ telescope to search for 
the 7=1-0, 7 = 2-1, and ./ = 3 - 2 rotational transitions of 
CF"^ toward the Orion Bar, a well-studied photodissociation re- 
gion (PDR) with an edge-on geometry favorable for the detec- 
tion of molecular species of small abundance. The observations 
are described in §2 below, and the results of those observations 
presented in §3. A discussion follows in §4. 

2. Observations 

The 7 = 1-0 rotational transition of CF"" at 102.58748 GHz 
(Plummer et al. 1986) was observed in two five-position strip 
maps in the Orion Bar, centered on the "Orion Bar (CO)" and 
"Orion Bar (HCN)" positions defined by Schilke et al. (2001; 
see our Figure 1 caption for coordinates). The 7 = 2-1 line 
at 205.17445 GHz and the 7 = 3 - 2 fine at 307.74438 GHz 
were observed at the central "Orion Bar (CO)" and "Orion Bar 
(HCN)" positions for each strip. 

The observations of CF+ 7=1-0 and 7 = 2-1 were car- 
ried out at the TRAM 30 m telescope, using the AlOO, BlOO, 
A230 and B230 receivers in position-switching mode, with 
an OFF-position located 10' away in azimuth. The backend 
employed was the VESPA correlation spectrometer; care was 

^ The Atacama Pathfinder Experiment (APEX) is a collaboration of 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observtory 



Fig.2. CF+ 7 = 1 - 0, 7 = 2 - 1 (IRAM 30 m) and 7 = 
3-2 (APEX 12m) spectra obtained toward the Orion Bar (CO) 
and Orion Bar (HCN) positions. (Spectra of lower signal-to- 
noise ratio, obtained in shorter integrations toward the Orion 
Bar (CO) position, were published previously by Neufeld et al. 
2005b.) 

Table 1. Observations of the Orion Bar 



Orion Bar (CO) Orion Bar (HCN) 

Line strength" ( J T^dv in mKkm s"') 
7= 1 -0 86 ± 10 72 ± 8 

7 = 2- 1 337 ±13 215 ±13 

7 = 3-2 428 ±34 183 ±26 

Line centroid (kms~' w.r.t. the LSR) 
7=1-0 10.5 ±0.15 9.6 ±0.15 

7 = 2-1 10.7 ± 0.05 9.8 ± 0.05 

7 = 3-2 10.6 ±0.1 10.0 ±0.1 

Line width (km s"' FWHM) 
7=1-0 2.7 ± 0.4 2.3 ± 0.3 

7 = 2-1 2.0 ±0.1 1.7 ±0.1 

7 = 3-2 3.0 ±0.3 1.5 ±0.2 

Column density (units of 10" cm"^) 
7=1 3.7 3.1 

7 = 2 6.0 3.8 

7 = 3 4.7 2.0 

Estimated total 19 11 



" Errors given are 1 sigma statistical errors. We estimate the sys- 
tematic uncertainties as ±20%. 

taken to move the known backend defects away from the ob- 
served lines. The CF"^ 7 = 3-2 transition was observed using 
the APEX 12m telescope, located at an altitude of 5100 m on 
the Chajnantor plateau in the Atacama desert of Chile, using the 
APEX-2a receiver. The observing mode was the same as at the 
30 m, with the same OFF-Position. As backend, the MPIfR Fast 
Fourier Spectrometer was used. Calibration at both telescopes 
is performed using a dual load scheme and an atmospheric 
model. The beam sizes were 24", 12" and 21" (HPBW), re- 
spectively. Since the emission is extended, in the following an- 
terma temperature units are used. 
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Fig. 3. Rotational diagram for CF"^ . The error bars include sys- 
tematic uncertainties, which we estimate as +20%. The solid 
curves show a fit to the data, in which the critical density for 
each transition is assumed proportional to the radiative decay 
rate. 

3. Results 

Figure 1 shows the CF^ 7=1-0 spectra observed at each 
of the positions in the Orion Bar, with the observed position 
for each observation indicated on a '^CO 7 = 3-2 map of 
the source. Emission in the 7=1-0 line is evident at 8 of 
the 10 observed positions. Like other optically-thin emission 
Hnes such as SiO 7 = 2 - 1 (Schilke et al. 2001), the CF+ 7 = 
1 -0 line peaks southeast (i.e. further from the Trapezium stars) 
than "CO 7 = 3-2. Rather than indicating that CF+ and SiO 
are located further away from the ionization front than '^CO, 
this behavior reflects the different sensitivity to temperature and 
density - and the different optical depths - of the lines from 
these three species. In particular, the '^CO 7 = 3-2 line is 
of moderate optical depth (t > 1), so its intensity probes the 
temperature at the t ~ 1 surface and is relatively insensitive 
to the total line-of-sight column density. Figure 2 shows the 
7 = 1-0, 7 = 2-1 and 7 = 3-2 spectra obtained at the 
central Orion Bar (CO) and Orion Bar (HCN) positions. Table 
1 lists the integrated antenna temperature, line centroid, and 
line width for each observed transition at these two positions. 

Figure 3 shows a rotational diagram for the two positions 
toward which all three transitions have been observed, com- 
puted for an assumed dipole moment of 1.07 D (Peterson et 
al. 1990). Here, we assumed that the line emission is opticaUy- 
thin, as suggested by the observed line strengths unless the cov- 
ering factor of the emission is implausibly small. The column 
densities for 7 = 1, 2 and 3 are given in Table 1. At both po- 
sitions, any reasonable extrapolation of the rotational diagram 
imphes that the 7 = 1,7 = 2, and 7 = 3 states are the three 
most highly-populated rotational states; thus the total column 
density we infer is only weakly dependent upon the extrapola- 
tion procedure adopted. The solid lines show the best fits to the 
rotational diagrams, which imply total CF^ column densities of 
A?(CF+) = 1.9 X 10'^ cm-2 and 1.1 x lO'^ cm'^ respectively for 
the lines-of-sight to the Orion Bar (CO) and Orion Bar (HCN) 
positions. Given our model predictions for the typical tempera- 
ture in the CF"^ emitting region, and as suggested by the curva- 
ture of the rotational diagrams, we assumed here that the slope 



is mainly controlled by the density rather than the tempera- 
ture, the former being insufficient for coUisional deexcitation 
of 7 = 3 to donunate spontaneous radiative decay. 

4. Discussion 

The results presented in Figure 1 and Table 1 indicate that CF^ 
column densities ~ 10^^ cm"^ are present over an extended re- 
gion within the Orion Bar. Such column densities are in good 
agreement with theoretical predictions (NWS), which predict 
CF"^ to be produced by the reaction sequence: 

Hz-hF^HF-hH (Rl) 

HF + C^ ^ CF+ + H (Rl) 
and destroyed primarily"* by dissociative recombination 

CF-' -I- e ^ C -I- F iR3) 

Because the reaction (Rl) is exothermic and moderately 
rapid even at low temperature (Zhu et al. 2002), HF becomes 
the dominant reservoir of fluorine in the gas phase, even rela- 
tively close to cloud surfaces. Beneath the surface of the UV- 
illuminated cloud, HF forms at precisely the point at which hy- 
drogen becomes molecular, and long before carbon gets incor- 
porated into CO (see NWS, Figure 2). Thus there is a substan- 
tial area of overlap between C"' and HF where reaction (R2) can 
produce CF'' abundances ~ few x 10"'"; here CF^^ accounts for 
~ 1 % of the solar abundance of fluorine. 

Figure 4a shows the CF'' column density predicted by the 
NWS model, as a function of the gas density, hh, and the in- 
cident radiation field xuv (again for the case of one-sided il- 
lumination at normal incidence). The results apply to the case 
where fluorine depletion is estimated using the simple treat- 
ment given in NWS §3.2 The predicted and observed results are 
in good agreement, although the remarkable agreement sug- 
gested by Figure 4a is undoubtedly fortuitous for several rea- 
sons: (1) Figure 4a compares the predicted column densities 
for a photodissociation region (PDR) viewed face-on with the 
observed column densities in a nearly edge-on PDR (in which 
optically-thin lines are significantly Umb-brightened); (2) the 
theoretical predictions are based upon several parameters that 
are substantially uncertain: e.g. the reaction rate coefficients for 
(R2) and (R3), and the gas-phase F abundance (see footnote 4). 

The CF'^ column densities plotted in Figure 4a are found to 
track the corresponding C column densities remarkably well, 
with the N(CF^)/N(C'^) column density ratio (Figure 4b) lying 
within a factor ~ 2 of 1.7 x 10"^ for every case represented in 
Figure 4. This behavior can be understood by considering the 
rates of CF^ formation and destruction via reactions (R2) and 
(R3). In equilibrium, we expect a ratio 

n(CF^) ^ ^2 «(HF) ^ J T \° ''«(HF)nc 

n(C+) ki n, \300K/ tip He 

Reaction with H2 to form HCF+, considered by Morino et al. 
(2000) to dominate the destruction of interstellar CP*, is in fact sub- 
stantially endothermic^ (by ~ 2 eV) and therefore of negligible impor- 
tance (NWS). 



4 



Neufeld et al.: Discovery of interstellar CF' 



00 

3 10.5 
10 

— -5.5 





r 10' 
: 10' 






; 10' 






: 10" 
-, , , 1 , , 




. 1 1 




' 10= 




. 1 . . . 


lO*"" -— 

1(? — 

1 


X 

1 1 . 



2 3 4 5 

Fig. 4. a) Upper panel: CF"^ column density, A^(CF"^), pre- 
dicted by the NWS model, as a function of the normalized 
incident radiation field Xuv- (b) Lower panel: corresponding 
N(CF^)/N(C*) ratio. Each curve is labelled with the assumed 
density, hh- The crosses represent the mean observed values 
for the Orion Bar (CO) and Orion Bar (HCN) positions, along 
with the radiation field = 3 x 10'* estimated by Herrmann 
et al. (1997) for the Orion Bar, while the open circles represent 
the corresponding model prediction. 



where k2 and kj are the rate coefficients given by NWS for (R2) 
and (R3), and nF/rin = 1.8 x 10 and ric/nn = 1.6 x 10 are 
the assumed gas-phase abundances of F and C. In the region 
where the and CF"^ abundances are large, ~ n{C^) ~ nc 
and n(HF) ~ np- 

Kuiper Airborne Observatory (KAO) observations of the 
C"^ ^^"3/2 Pi/2 (158//m) fine-structure line carried out toward 
the Orion Bar (Herrmann et al. 1997) permit us to compare our 
theoretical prediction for A^(CF^)/A^(C^) with the measured 
value. Since both the C"^ and CF^ transitions are optically- 
thin, such a comparison has the merit of eliminating any hmb- 
brightening effects, all the observed lines being enhanced by 
the same factor. Based upon their observations of the Orion 
Bar, Herrmaim et al. inferred a Une-of-sight C"^ colunm density 
of 3 X 10'^ cm"^, averaged over a 55" beam. Given our esti- 
mates of A^(CF^) at two separate positions along the midline of 
the bar, we may estimate the N(CF'^)/NiC'^) ratio as ~ 5 x 10"^, 
a factor of ~ 4 smaller than the theoretical prediction. This dis- 
crepancy may suggest (weakly, in light of all the uncertainties) 
that the F abundance assumed by NWS is too large and/or that 
the reaction (R2) proceeds less rapidly than was assumed. 

Despite these uncertainties, our detection of CF"^ lends sup- 
port to our theoretical model of fluorine chemistry. In partic- 
ular, the observed CF"^ column density argues strongly for a 
substantial region of overlap between and HF. This, in turn, 
supports our contention that HF becomes abundant very close 
to cloud surfaces, and suggests that HF may prove a valuable 
tracer of molecular material in difl'use clouds. Observations of 
HF and CF"^ may also provide an important tool for studying 
the elemental abundance and nucleosynthetic origin of fluo- 



rine^. As noted by NWS, absorption line observations of the 
1232 GHz UF J = 1-0 transition will be possible with the 
Herschel Space Observatory and the Stratospheric Observatory 
for Infrared Astronomy; in addition, observations of HF at cer- 
tain redshift ranges beyond z = 0.3 are possible at ground- 
based observatories and will be greatly facilitated by ALMA. 
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' The elemental abundance of fluorine is of particular interest be- 
cause the nucleosynthetic origin of fluorine remains controversial. 
Fluorine is produced in stars on the asymptotic giant branch (Jorissen, 
Smith, & Lambert 1992; Werner, Rauch, & Kruk 2005), but there is no 
consensus on the relative contributions of other sources such as Wolf- 
Rayet stars and type 11 supernovae (e.g. Renda et al. 2004; Palacios,, 
Arnould, & Meynet 2005). Because AGB stars exhibit a large en- 
hancement in the fluorine abundance (up to a factor ~ 250), it will 
be interesting to search for CF* in circumstellar envelopes as well as 
in the wind-blown bubbles surrounding Wolf-Rayet stars (Marston et 
al. 1999, Marston 2001, Rizzo et al. 2003). 



